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This paper and the results therein are geared towards building a basic toolbox for 
calculations in quantum information theory of quasi-free fermionic systems. Various 
entropy and relative entropy measures are discussed and the calculation of these 
reduced to evaluating functions on the one-particle component of quasi-free states. 

The set of quasi-free affine maps on the state space is determined and fully charac- 
terized in terms of operations on one-particle subspaces. For a subclass of trace pre- 
serving completely positive maps and for their duals, Choi matrices and Jamiolkowski 
states are discussed. 
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1. INTRODUCTION 

There are not too many classes of states or quantum operations that can be handled in 
detail. Well-known examples are gaussian structures in bosonic systems [1] and exchangeable 
states and maps for spin systems [2, 5]. Although the systems under consideration can often 
be quite large, even infinite, the computational difficulty of most associated quantities is 
many orders of magnitude lower than what one encounters in more general systems. In 
particular gaussian states have been used in quantum optics to that effect for many years. 
Perhaps because of the close link between the two fields, gaussian states were also the first 
of the aforementioned classes to come under consideration in quantum information theory. 
Recent years have seen a large amount of work done on their role in bosonic systems, see 
for instance [7]. 

This paper deals with quasi- free fermionic systems [1]. In field theory and statistical 
mechanics such effective free evolutions and states have been used extensively as an approx- 
imation to interacting systems, a well-known example being the Hartree-Fock approximation. 
The main simplifying feature lies in the particular combinatorial properties of correlation 
functions and maps. In fact, states and maps are fully determined by one-particle opera- 
tors. As the dimension of the observables increases exponentially with the dimension of the 
one-particle space we obtain a very significant reduction of complexity. 

Fermionic systems and quasi-free states should be of particular interest to information 



* Electronic address: 
^Electronic address: 

* Electronic address: 



brccht.dicrckx@fys.kuleuven.be 
mark . f anncs@fys .kulcuvcn . be 
dokmpo@univ.gda.pl 



2 



theorists. A qubit system can always be mapped on an interacting system of fermions by 
the so-called Jordan- Wigner isomorphism. A particular subset of quantum operations 
on qubits can then be identified with quasi-free evolutions of fermionic systems and in these 
problems quasi-free states play an important role. A recent article dealing with this duality 
between fermions and qubits is [6]. 

We attempted to write a rather self-contained paper, offering a toolbox for computations 
and testing ideas in quantum information. Exponential elements instead of the standard 
creation and annihilation operators, see [3], appear to be an efficient computational tool. 
Well-known objects, as quasi-free states, are reconsidered in these terms but the main goal 
are quasi-free quantum operations. As we restrict ourselves on purpose to finite dimensions 
we only need linear algebra. This limitation can be overcome: infinite fermionic systems no 
longer have a canonical representation and so one has to use appropriate representations. 
This typically involves introducing trace-class conditions on the one-particle operators. 

The paper is organized as follows: section 2 recalls some basic properties of fermionic 
Fock space and creation and annihilation operators satisfying CAR (canonical anticommu- 
tation relations). In section 3 we introduce the basic exponential operators and study their 
properties. Section 4 reconsiders quasi-free states and introduces new calculation techniques 
for some entropic quantities. Section 5 deals with quasi-free quantum operations. 



2. FERMIONIC FOCK SPACE 

The quantum mechanical description of fermions is firmly connected to the mathematical 
concept of antisymmetric Fock space, especially in finite dimensions. Throughout the article 
we will regularly fall back on this to prove, calculate or physically motivate our expressions. 

In the following, will denote a finite dimensional complex inner product space, which to 
almost all intent and purpose, we can regard as a finite Hilbert space. Many of the results 
can be extended to infinite dimensional Hilbert spaces, modulo, of course, some suitable 
additional conditions. 

The symbol ® is well-known as a notation for the tensor product of two or more objects. 
In this, we will use the symbol A to denote the antisymmetric tensor product of vectors, 
operators and even algebras. Although the meaning of the symbol will change depending 
on the setting in which it is used, restricting ourselves to a single symbol, greatly simplifies 
the notation and looks, at least to a physicist's eye, more elegant. As is often the case with 
degenerate notations, the context should specify which version of the wedge we are talking 
about. 

2.1. Antisymmetric Vector Spaces 

We define the fc-fold antisymmetric tensor product or wedge product of vectors 
(ft, ip 2 , ■ ■ ■ , ¥k in F) as 



(fx A y? 2 A . . . A ip k ■= 




e(cr) (/^(i) <g> </V( 2 ) <E> • • • ® iPa(k) 



where a runs over all permutations of the k indices and e(a) = ± depending on the parity 
of the permutation, + if even, — if odd. 
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Let U a be the unitary operator which implements the permutation a on (£) k $j. The fully 
antisymmetric subspace of (3 k $) consists of these vectors 7] G (E) k $j~ which satisfy U a i] = e(a)r]. 
It is spanned by the A;-fold antisymmetric vectors and we denote it as A fc 5) or $)( k \ This 
space has the same inner product as which can be written more succinctly as 

(p 1 A p 2 A ••• Ap k , A^ 2 A ••• A^ fc ) = Det ([(</?;, ^j)]^-)- 

This in turn makes it quite easy to transport a basis {ei, e 2 , . . .} of S) to To each such 
orthonormal basis corresponds an orthonormal basis 

{e A : Ac{l,2,...}&#(A) = fc} (1) 

of fj^ where 



e A := e h Ae i2 A ■ ■ ■ Ae ik , i x < % 2 < ■ ■ ■ < ik & A = {i u % 2 , . . . , it}- 

The wedge product extends naturally to antisymmetric vector spaces. If (p e i^ ( - fc - ) and 
then 

^Ax:= ; V e(A, M)ip A ® X m 



is an element of ^ h+t \ The sum runs over the ordered partitions of {1,2, . . . , k + £} in 
subsets A and M with #(A) = k and #(M) = £, v?a and xm are the ordered injections of 
ip and x in the corresponding tensor products and e(A, M) is the parity of the permutation 
(A, M). To clarify this a bit, assume that tp and x are elementary antisymmetric tensors, 
i.e. 

p = pi A . . . A p k and x = Xi A • • • A Xe- 
Then there is a canonical way to define the wedge product between these two vectors as 

<P A X = <Pi A • • • A p k A Xi A . . . A Xi 

= in tk 1 , Ve(ff)[/ g (^i® ...®<Pk®Xi® ■■■Xe)- 

The above defined extension to wedge products of general antisymmetric vectors is then just 
the linear extension of this and so it becomes easy to prove the following important property 
of the wedge. 

Lemma 1. The wedge operation is associative 



(p A X ) A^ = p A ( X Aij). 



2.2. Construction of the Fermionic Fock Space 



Although the A;-antisymmetric vectors we introduced in the previous section are an apt 
description for the states or wave functions of k fermionic particles, we need more to properly 
describe a physical system of such particles. In nature there are processes which do not 
conserve the number of particles. So we need a setting in which we can jump between 
different A;-antisymmetric vector spaces. 
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Consider for example a fermion with d modes. The state space associated with it is 
then also ^-dimensional. If we add another identical particle to our system, the degrees 
of freedom do not go up as they do with qubits, rather, since we need to obey the Pauli 
exclusion principle, the dimension of the state space, which is the 2-antisymmetric space, 
is only ( ). In general, the dimension of the state space of k such particles is as can 
be readily seen from (1), and in particular, the state space of d such particles is only one- 
dimensional. So we can only combine d fermions before we literally run out of space to put 
them. 

To unify these concepts of particle creation/annihilation and the exclusion principle, the 
antisymmetric Fock space is introduced. The fermionic Fock space with one-particle space 
Sj is 

T{S}) := C©-f>©f> (2) ©•••©€. 

It follows from the paragraph above that T(7i) has dimension 2 d where d = dim(j}). The 
first term in the direct sum is the vacuum state and the last one is the completely filled 
Fermi sea. 

Piecing two systems together 

Given an orthogonal decomposition of $) into subspaces Sji and fj 2 , there is a natural 
isomorphism between and the tensor product of the Fock spaces with fji and Sj2 as 

one-particle spaces 

r(£i©£ 2 ) =r(#i)®r(J3 2 ) 

explicitly given by 

(</?i © Vi) A • • • A (if k © i) k ) = A • • • A ip ir ) © {i) n A • • ■ A i) ja ). 

The summation sign J^® points at mixed sums and direct sums and the summation runs over 
all ordered partitions of {1, 2, ... , k} in two subsets {zi, 12, ■ ■ ■ , v} an d {ji, J2, • • • ,js}- Due 
to the antisymmetry one has to pay attention to the order of the factors in this isomorphism. 

Elementary antisymmetric vectors 

We would also like to point out a peculiar property of elementary vectors which will come 
in handy later on. 

Lemma 2. A nonzero vector ip E is an elementary vector, i.e. can be written as 

(p = ijj x A . . . A Vfe 

if and only if the space 

{ x ef)\ X A V = 0} 

is k- dimensional. 

Proof. Consider an elementary vector 95 in Sj^ which can be written in the simple form 



if = ipi A . . . A ip k 



5 



where {ipi}i U {xe}e forms a basis for the generating space 5} and {ipi}i^{xe}i- From the 
above, it is clear that a set of vectors is linearly dependent if and only if the elementary 
tensor constructed from them is zero. The xi are linearly independent of the ipi and so 
cannot contribute to the set in the lemma. Only vectors which are built up exclusively out 
of the ipi contribute and the space generated by the ipi is ^-dimensional. 

Now assume that the set mentioned in the lemma is indeed dimensional. This implies 
that we can find (d — k) orthogonal vectors in the ci-dimensional 5} which are linearly inde- 
pendent of this set and thus of the constituting vectors of i p>. So only k linearly independent 
vectors can be involved in the creation of (p. But the only non-zero fc-antisymmetric vectors 
which can be built out of k vectors, are proportional to each other and elementary. □ 

2.3. The CAR Algebra 

Fermionic Fock space can also be built in another way. Any vector (p of 53 induces a 
linear operator a* ((f) from 53^ to fj( fc+1 ). We first define the action of a*(cp) on elementary 
vectors and then linearly extend it to the whole space. 

a*0)(V>i A ■ ■ ■ A Vfc) := <P A V>i A . . . A ip k . 

The operator a* (if) is called a creation operator, its adjoint an annihilation operator and in 
effect they emulate the creation or destruction of a Fermion in the state p>. By repeatedly 
applying creation operators to the vacuum vector, we can build up the entire Fock space. 

Although we will not use this language very frequently, the operators defined above satisfy 
exactly the CAR required for a quantum mechanical description of Fermions 

{a(^),a(^)} = and {a(cp), a*(V>)} = (0 , V) L (2) 

The algebra built on the creation and annihilation operators is called the CAR-algebra 
in reference to the important commutation relations (2). It coincides with the algebra of 
linear transformations of the Fock space r(53). This algebra is in fact a universal algebra 
because it is the unique algebra generated by a unit element 1 and by {a(<p) : <p G 53} such 
that the operators a(ip) satisfy the CAR-conditions and that the map ip \— > a(ip) is complex 
antilinear. We will denote it by 21(53). 

The following proposition will be needed later on. It is well-known, so we state it without 
proof, see [1]. 

Proposition 1. Suppose that uo is an even state, i.e. vanishes on monomials in creation 
and annihilation operators with odd number of factors, on and that a is a state on 

then there exists a unique state uj A a on 21(53) A 2l(£) := 21(5} © defined by 

(u A a) (xy) := u){x) <x(y), x G 21(5}), y G 2l(£). 

Remark 1. The wedge product 21(53) A 2l(£) we implicitly defined in the above proposition 
is not the same as the tensor product of the two algebras as can be seen from the following 
constructive explanation. 

There is a natural embedding 

3l : 2l(53i) ^ 2l(53i © 53 2 ) : 3i{a(<p x )) = a{p 1 © 0), Vl G 53i 
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and of course an analogous embedding j 2 of 2l(fj 2 ). Clearly, jx($l($)i)) and 72(21(^2)) gen- 
erate 2t(f)i © ^2) but they do not sit in © fj 2 ) as tensor factors because 

{a # (<^i©0), a # (0ffiy? 2 )} = instead of [a # (^ © 0) , a # (0 © y? 2 )] = 0, 

a # denotes either a or a*. Sometimes, © i} 2 ) is called the graded tensor product of 

Sl(i5i) and 2t(ij 2 ). 



3. THE GICAR ALGEBRA 

For any one-particle basis {ej}, a special operator can be defined 

X:=^a*(e>(e,). (3) 

i 

It is invariant under the gauge group of 21(f)), but also under any arbitrary basis transforma- 
tion of the one-particle space and so (3) defines a unique operator in the algebra. Consider 
the action of N on an arbitrary (non-zero) fc-particle vector 

N (fx A . . . A ifk = k (fx A . . . A ifk- 

So, iV counts the number of particles in a given state and as such is called the number 
operator. Its eigenspaces are obviously the fc-antisymmetric spaces and its spectrum consists 
of the integers {0, . . . , d}. The commutant of the number operator, is called the gauge 
invariant CAR-algebra or GICAR for short. It is the largest subalgebra of 21(f)) that is 
invariant under the gauge group. It is also generated as the span of all monomials in a,a* 
containing as many a's as a*'s. 



Exponential elements 

With d = dim(fj) and the algebra of complex square matrices of dimension k, the 
GICAR-algebra can be written as 

2l GICAR (fj) = C © M d © ©•••©€. 

The dimension of 2l GICAR (fj), counted as a complex vector space, is then seen to be ( 2 ^). 

Remark that the GICAR is the subalgebra of the CAR of block diagonal transformations 
of Fock space, in particular it contains elements of the form 

E(A) := 1 © X © (X © X)\ m © • • ■ © (® d X)\ m 

or with an obvious notational meaning 

E(A) = 1©A©(AAA)©---© (A d X). 

Their spectrum <t(E(X)) can be computed to be 

(j(E(A)) = \Y[\ : A; G ct(X)&A C {l,...,Rank(X)}| 
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This property is easily verified by looking at the antisymmetric part of ®fcX. As a tensor 
product, the eigenvalues of this are exactly all possible monomials of the eigenvalues of X 
of length k. If we number the eigenvalues of X repeated according to their multiplicities as 
Aj, then a monomial of the Aj will contribute to the spectrum of A^X if and only if each Aj 
appears exactly once or not at all. 

They also enjoy the following properties 

E(l) = 1 

E(xy = Epr) 

E(X)E(Y) = E(XY) 
E(X) > iff X > 
E(X 1 ®X 2 ) = E(X X ) ® E(X 2 ) 
TrE(X) = Det(l + X) 

2l GICAR (£)=Span({E(X)}). (4) 

The first six properties can be easily checked by looking at the spectrum of the E-operators 
and the isomorphism we mentioned above. The last property follows from expanding 
A i — > E(AX) around A = and remarking that Span({(g> fc X|^ {fc) }) coincides with the linear 

transformations of $j( k \ 



k-P article projectors 

Another useful set of elements consists of certain one-dimensional projectors. Given a 
/c-dimensional subspace & of fj, all vectors <pi A • • • A fk with (fx, . . . , iff. G & are proportional 
and span therefore a one- dimensional subspace of $)( k \ The projector on that space will be 
denoted by P*(&). If we denote by [M\ the projector on A then 

P*(£) = [£] ® [&} w . 

It is convenient to associate the projector on the vacuum space with P*(0) where is the 
zero-dimensional vector space. 

These projectors are contained in the closure of Span({E(X)}) and as such inherit all 
relevant properties of the E(X) listed in (4). They arise as the limits of normalized E- 
operators, 




< X n < 1. 



Remark 2. The map E defined by 

E(X) ^ Det(l -X) E (j^^)> V0<X<1 

is uniformly continuous and extends therefore continuously to [0, 1]. This extension sheds 
some light on how the n-particle projectors arise as limits of E-operators. 

Suppose that 1 is a fc-degenerate eigenvalue of X. X can then be decomposed as a direct 
sum of a projector P and a (d — /c)-dimensional object: 



X = P@X 
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;i-6„) fc Det(l-X) E(-i=-P 



For any sequence < (e n ) n < 1 that converges to one, the operators 

E(e n P©X) 

are well-defined, bounded and as per (4) isomorphic to 

l-e„ / VI -X 
By rearranging the factors in this expression, we get 

E(e n P)®E(X). 

We can write out the first factor as 

(1 - e n ) k { (0(1 - e n )^V n a' P) (0 0) }, 

and then it is clear that this will converge to 

© ... © (A k P) © ... © 0. 

We will often ignore the possibility that 1 is included in the spectrum of X when we 
are talking about expressions containing E(X). When 1 is contained in the spectrum of X, 
/(E(X)J should be interpreted as 

lim f(E(X n )), 0<X n <l 

x n — >x 

whenever / is a continuous function. 

Remark 3. We can express a general E(X) using the above projectors as 

E(X) = ^x A P*(£ A ) 

A 

where A plays the same role as in (4) and the x\ are the products of the corresponding 
eigenvalues. So in essence, this gives us the eigendecomposition of E(X). 



4. QUASI-FREE STATES 

The chief objects under study in this article are the so-called quasi-free states and maps. 
These states are the fermionic counterparts of gaussian measures for classical systems or 
gaussian states for bosonic systems. They are sometimes referred to as determinantal states 
or processes for reasons that will soon become obvious. 

A linear functional u on the CAR algebra which assigns zero values to all monomials in 
creation and annihilation operators except for 

u(a*((pi) ■ ■■a*{(p k )a{ip k ) ■ ••aW) = Det f [(^ , Q(fj)]j 

extends to a state on 21(f)) if and only if the linear one-particle space transformation Q 
satisfies < Q < 1. Such an u is called gauge invariant quasi-free and Q its corresponding 
symbol. The notation ujq will be used to connect the state to its symbol. 

Using the language developed in section 3 we can calculate the density matrix pg corre- 
sponding to a state ujq. 
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Lemma 3. The density matrix pq corresponding to a state ujq with symbol Q can be written 
down explicitly as 

E(Q) = Det (1 - Q) (l © -Q— © ( A Q 



1-Q \1-Q 1-Q 

Proof. Consider the symbol Q of a general quasi-free state to. We can always find a ID 
decomposition of Sj such that it is amenable with the eigendecomposition of Q, i.e. 

Q = y^gi|ei)(ej[) = C |ej) & ©» = i^- 

j 

By straightforward computation we can check that cuq = A^o^ with Qi : = ^|e»)(ei| and 
since 

r(£i © i3 2 © • • • © sj n ) = r(£i) © r(# 2 ) © • • • © r(J5 n ), 

= ®iSt(«i). 

For the one- dimensional Hilbert space fjj, the Fock space r(f)j) is 2-dimensional and the 
associated representation can be expressed on C 2 as 

M°) = (i o) ' \ Qf ^ = (i 
The relevant part of Q on this Fock space is then 

Qi = (o l-q 

Because of Proposition 1 and the uniqueness implied in there, the wedge state Aj ojq i 
must be isomorphic to the product state 



We can rewrite this in terms of exponential elements by 

1 - qij \0 l-q d 

Det ° jj, ft^l 



Det (1 - Q) E 



Ql \ _ „ r / 



which is isomorphic to 



p Q = Det (1-Q) E 



Q 



1-Q 



Finally, by continuity of the map E this results also holds for = 1 and so in particular for 
projectors. □ 



10 



By remark 3, there is an alternative way to write this density matrix. Let < Q < 1 be a 
linear transformation of the one-particle space $) and let {ei, e%, . . . , e^} be an orthonormal 
set of eigenvectors of Q, i.e. 

Q e i = Qj e i> <qj <1, j = l,2,...,d. 
For a subset A of {1, 2, ... , d} define 

q A : = JJg r ] j (1 - q s ) and fi A := Span({e r : r G A}) 

reA se{l,...,d}\A 

then 

= qA P * (^ A ) ' 

A 

As the P*(S)\) project on mutually orthogonal subspaces, this decomposition shows that the 
q A are the eigenvalues of pq. This allows to explicitly compute quantities as the Renyi and 
von Neumann entropies of Pq. 

As every P*(^) defines a pure state on 2t(i3) that is gauge invariant and quasi-free, we see 
that every gauge invariant quasi-free state is a convex mixture of pure gauge invariant quasi- 
free states. The projectors P* are generally a small subset of the one-dimensional projectors 
acting; on This can be seen by a simple parameter count. In order to parameterize a 
generic m-dimensional complex subspace of C n , we need 2m(n — m) real parameters. As 
dim(i)W) = g) with d = dim(£) we need 2 ((f) — l) real parameters to specify a generic 
one-dimensional subspace of that is to say a /c-particle pure state on 2l GICAR (fj) while 
we need only 2k(d — k) < 2((*?) — l) real parameters to specify a fc-dimensional subspace 
of fj which corresponds to a pure quasi-free state. Therefore, the convex hull of the gauge 
invariant quasi-free states is strictly smaller than the state space of -4. GICAR (5}). The linear 
span of the pure quasi- free states coincides however with all linear functionals on A GlcAR (Sj). 

To illustrate that the quasi-free states do not form a convex set on their own, we provide 
the following proposition which will also be of use later on. 

Proposition 2. Let ujq 1 and ujq 2 be quasi-free and letQ<X<l, then \ojq 1 + (1 — A) uq 2 
is quasi-free iff Qi — Q2 is of rank or 1. Moreover, if the rank condition holds, 

^ u Qi + (1 - ^) U Q2 = ^AQ!+(l-A)Q 2 - 

For the proof of this we refer to [8]. 



4.1. Entropy Related Quantities of Quasi-free States 

The knowledge we have about the eigenvalue decomposition of a quasi-free state ujq 
allows us to translate the general expressions for entropy related quantities of that state 
to expressions on the one-particle density matrix. Especially calculations involving p- Renyi 
entropies become very simple in this way. 

Proposition 3. The p-Renyi entropy of a quasi-free state ujq is 

H p (uj q ) = -L_ Tr log((l - Qf + Q p ). (5) 
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Proof. 



1 


1 




p 




1 




1 




p 




1 




1 




p 




1 




1 




p 



log |Det + (jzTqJ ) Det (1 - Qf 
logDet((l -Q) p + Q p ) 
Trlog((l-Q) p + Q p ). 



Proposition 4. The von Neumann entropy of a quasi-free state uoq is 

S(w Q ) = -Tr(QlogQ + (l-Q)log(l-Q)' 



□ 



Proof. The von Neumann entropy is the limit of the p-Renyi entropy when p I 1. As the 
expression in (5) becomes indeterminate (0/0) we use de l'Hopital's rule combined with 
Jacobi's formula for differentiating a determinant 



This yields 



^Det (A) = Tr(A- ld ^) Det (A), 
ax \ ax J 



S(u Q ) = -lim^- log Det ((1 - Qf + Q p ) 
p|i dp 

= -lim^-Det((l-Q) p + Q p ) 

pU dp ' 

= - Tr (Q log Q + (1 - Q) log(l - Q)) 



□ 



Proposition 5. Let Q\ and Q2 be two symbols such that ker Q2 C ker Q\ and ker(l —Q2) C 
ker(l — Qi), then the relative entropy of ujq 2 with respect to ujq 1 is given by 



s (u Qi ;uq 2 ) = Tr|Qi(logQi -logQ 2 ) 

+ (1 - g0 (log(l - Q x ) - log(l - Q 2 )) } . 

Proof. As above, the computation uses an appropriate limit 

S(p; a) := Tr p(logp — loger) = lim ^— Trf// +1 — pa p \ 

p^o dp V / 
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Using (5) 

= lim ^ {Det (Q? +1 (l - Q i y+ 1 ) - Det ((1 - Q x ) (1 - Q 2 f + Q l Q p 2 ) } 
= TrjQi (log Qi - log Q 2 ) + (1 - Q x ) (log(l - Qi) - log(l - Q 2 j) }. 

□ 

5. QUASI-FREE COMPLETELY POSITIVE MAPS 

Let A be a completely positive map on Aid- Dual to A is another completely positive 
map A* on Aid 

Tr A(er) X = Tr a A*(X), a,XeM d . 

A general quantum operation may be described either in Schrodinger or in Heisenberg pic- 
ture. In the first case we use completely positive maps A with the additional property that 
Tr A(ct) = Trcr. Such maps restrict to affine transformations of the state space of Aid and 
are called trace-preserving completely positive (TPCP). Their duals A* leave the identity 
untouched and are therefore called unity-preserving completely positive (UPCP). 

5.1. Unital Maps 

We will consider here two families of UPCP maps on 21(5}) which generalize the expres- 
sions for quasi-free states [4]. Any unitary U on one-particle space defines an automorphism 
of through 

a* (ip) I— > a*(U ip). 

One can either check that the CAR are preserved or explicitly compute for an m-particle 
vector t/> 

E(U)a*{<p)E(U)*if) = E(U)a*(<p) ( U* <8> • • • g> U\ ip) 

m times 

= E(u)<ph (u* ®-..®tr^) 

= E(U)(U* ® ... ® [/*)(([/ v?) A V) 
= (U (p) Aip = a*(U<p)ip. 

Such automorphisms are called quasi-free. 

Next consider a couple A, of linear transformations of Sj such that 



< B < 1 - A* A. 



(6) 
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The block matrix 

A VI - AA> 



V :-- 



VI - A* A A* 



is unitary on fj © fj and, using the constraint on (A, B), we can always find a linear map 
< Q < 1 on fj such that B = a/1 — A*AQy/l — A* A. Using these ingredients, we 
construct a UPCP map A* AB on 21(5}) by concatenating three UPCP maps: the injection 
a* ((f) i— > a*(<y9 © 0), the quasi- free automorphism defined by the block unitary V and the 
projection idAujQ from 21(f) © fj) to 21(f))- An explicit computation yields 

a Ab( q *(^i) ' ' ' a *(fk)a(ipe) ■ ■ ■ a(ipi)) 

= y^ea> g (a*((p mi ) ■ ■ ■ a*(y? m Ja(?/v) ■ ■ ■ a(ip ni )j 

x a*{Aip h ) ■ • • a*(Aip ik _ r )a(Aijj jt _ r ) ■ -a(Aij) h ). 

Here the summation is taken over all ordered partitions 

{{ii, . . . ,ik- r }, {mi, - - - ,m r }} and {{ji, . . . ,je- r }, {ni, . . . , n r }} 

of {1, . . . , k} and e is the parity of the corresponding permutation. It is not hard to show 
that the condition (6) is also necessary to have B UPCP. This follows already from 
the requirement that any quasi-free state on 21(f)) should be mapped into a state, which 
incidentally will also be quasi-free. From the definition one sees that the map A* restricts 
to a UPCP map of ^l GICAR . 

To define the second family, we first need a complex conjugation: fix an orthonormal 
basis {ei, e%, . . . , e^} in C d . The elements of the basis will be considered as real vectors and 
the complex conjugation is the conjugate linear operator 

j j 

From this definition we see that (jp) = <p and ip) = (ip , <p). We also introduce the 
conjugate of a complex linear transformation A by Acp :— Alp. The transformation A is 
complex linear and satisfies 



A + aB = A + aB, (A) = A, A* = (A)*, and AB = AB. 
The entries of A in the distinguished basis {ei, e^, ■ ■ ■ , e^} are 



(A)ij — (ei , Aej) — (ei , Aej) — (Aej , e,) — (e^ , Aej) — A^ 



In particular the conjugate coincides with the transpose for Hermitian elements. 
The second family of maps we will consider is of the form 

r* A , B { a *M ■ ■ ■ a*(ip k )a(ijj e ) ■ ■ • a(V>i)) 

= ^euj B (a*(<p mi ) ■ ■ ■a*(<p mr )a('ip nr ) ■ ■ ■ a(ip ni )j 

x a(ATp-) ■ ■ ■ a(AJpi i r;)a*(A'ip je _ r ) ■ ■■a*{A^ jl ) 
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with the same summation convention as above, T* A B is UPCP if and only if < B < 
1 — A T (A T )* . 

Both Aa,b and Ta.b map gauge-invariant quasi-free states into gauge-invariant quasi-free 

states 

^a,b{uq) = uja*qa+b and T a , b {ujq) = U-a t q t (a t )*+b+a t (a t )*- 
We now compute the actions of B and T* A B on elements of the type E. 

Lemma 4. Suppose that < B < 1 — A* A and 1 — B + XB invertible, then 

A* AB (E(X)) = Det (1 - B + XB) 

xE(l + i(l-B + XB)-\X - 1)A* 
= Det (1 - B + BX) 

x Efl + A(X - 1)(1 -B + BX)- 1 A* 

Proof. Suppose first that 0<B<t — A*A, the general case follows by continuity. Choose 
now < Q < 1 and put Q' := ^(Q) = A*QA + B, then also < Q' < 1. Recall also that 
the density matrix pg = Det (1 — Q) E{Q/{1 — Q)) (and the analogous expression for pqi). 
We now compute using 

Det (1 + CD) = Det (1 + DC), (7) 
Tr(A AB (p Q )E(X)) 

( A*OA -\- B 
= Det (1 - A*QA - B) Det 1 + ^ -X 



1 - A*QA - B' 
= Det (l - A*QA -B + (A*QA + B)X) 
= Det (1 - B + BX + A*QA(X - 1)) 

= Det (1 - B + BX) Det (l + A*QA(X - 1)(1 -B + BX)' 1 ) 
= Det (1 - B + BX) Det (l + QA(X - 1)(1 - + BX)" 1 ^*) 
= Det (1 - B + BX) Det (1 - Q) 

= Det (1 + B(X - 1)) Det (1 - Q) 

x Det (l + — ^— (1 + A(X - -B + BX)~ l A*) 
= Tr(p Q A^(E(X))). 
The second form of the expression follows from (7) and 

(1 - B + XB)-\X - 1) = (X - 1)(1 - B + BX) -1 . 



□ 



Although the previous theorem completely defines the map A*, it might be useful to state 
how the action of this map looks on the density matrix of a quasi-free state. 
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Proposition 6. Let Q be the symbol of a quasi-free state ojq and suppose that 1 — Q + {2Q — 
1)B is invertible, then 

AX B (Po) = Det(l-g + (2Q-l)B) 

x E(l + A(t-Q + (2Q - 1)B)-\2Q - 1)A*). 

Proof. It suffices to replace X with in the proof of the previous proposition and add 
the required normalization. The result also remains valid when Q becomes a projector, or 
even more generally, when 1 G cr(Q) by the continuity of the map E. □ 

The maps T* can be handled in a similar way. 

Lemma 5. Suppose that < B < 1 - A T (A T )* and that 1 - B T - AA* + X T B T + X T AA* 
is invertible, then 

T\ B (E(X)) = Det (l - B T - AA* + B T X T + AA*X T ) 

x E(l + A(l - B T - AA* + X T B T + X T AA*)-\l - X T )A*^j . 

Proposition 7. Let Q be the symbol of a quasi-free state uoq, < B < 1 — A T (A T )* and 
suppose that 1 - Q T + (2Q T - 1)(AA* + B T ) is invertible, then 

^*a,b(pq) = Det (1 - Q T + (2Q T - 1)(AA* + B T )) 

x E(l + A(l - Q T + (2Q T - 1)(AA* + B T )y\l - 2Q T )A*). 

5.2. Trace-preserving Maps 

The stability of the set of quasi-free states with respect to quasi-free CPTP maps can 
essentially be used as a characterization of such maps. 

Proposition 8. The set of quasi-free states is invariant under a linear map T if and only 
if the action of the map T can be expressed as 

r (Po) = Pi(Q) (8) 

where j(Q) = ±A*QA + B or j(Q) = ±A*Q T A + B. Furthermore -f(Q) = A*QA + B is 
CP iffO<B < 1 - A* A and i(Q) = —A*Q T A + B is CP iff A* A < B < 1. 

Proof. Consider a map 1 T : (% GICAR )* _> (yGiCARy with the f n ow j n g properties: 

• trace-preserving 

• C-linear 

• maps quasi-free states onto quasi-free states 



The * in this formula denotes the dual of the algebra, i.e. the span of the state space. 
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Since the map T leaves the set of quasi-free states invariant, there must be a corresponding 
map on the set of symbols, so we propose a map 7 as in (8). Consider now a matrix Q which 
is an element of the open interval [0, 1[ in 9Jl(C d ) and a one-dimensional projector P on the 
same algebra. If e is small enough (but non-zero) , Q + eP is still an element of the unit 
interval. Since the difference of Q and Q + eP is of Rank 1 

^Q+XeP = ^(l-X)Q+X(Q+eP) = (1 ~ A)wq + \UQ +tP . 

And because of the linearity of T, this gives us the following results about 7 

Rank { 7 (Q + eP) - 7(g)} = or 1, Ve 

and 

(1 - A)t(Q) + A 7 (Q + eP) = 7 (Q + AeP). 
We can rewrite this in a more useful form 

l(Q + eP) = ^{ 7 (Q + AeP) - (1 - A) 7 (Q)} 

= ef {7(Q + AeP) -7(g)} + 7(Q) 
Ae 

= 7 (Q) + ed(Q,P) 

where in the final line we have introduced the function d(Qi, Q2) with the property that if 
Q2 is a one-dimensional projector, the Rank of d is or 1. 

As the maps T, E and E^ 1 are all Frechet differentiable on the unit interval, 7 should also 
be differentiable. The function d we introduced above, is actually the derivative of 7 and 
as such it is linear in the second argument. So clearly for any matrix < Q = g^Pj < 1 
such that the Pj are one- dimensional projectors, 

t(q) = * p <) = J2 * p *) + ^^°) 

i j 

Because of the Rank conditions on d(Q, P) condition, we get 

7(Q) = B± A*QA or 7 (Q) = P ± A*g T A. 

The conditions < P < 1 - AM for the case 7(g) = A*gA + P and A* A < B < 1 for 
the case 7(g) = — A*Q T A + B are necessary because states have to be mapped onto states, 
in particular the image of a symbol should remain a symbol. Conversely, if the conditions 
hold, then we may define the dual UPCP maps as in section 5.1 □ 



5.3. Jamiolkowski states and Choi matrices of quasi-free completely positive maps 



There are two ways of encoding a CP map, in fact of encoding any super-operator, going 
under the names of Jamiolkowski state and Choi matrix. Let {e^} be the standard basis of 
C d with associated matrix units := |ej)(ej|. The Jamiolkowski state J(r) of T is defined 

as 



J(r) := (id®r) 



Z^ 1 



E 



7 , e *i 



r(e, ; 



17 



It is easily verified that J(T) is a state whenever T is TPCP, i.e. a quantum operation in 
Schrodinger picture. The Choi matrix of a CP map T is a very similar object 

mostly used in Heisenberg picture. So if T is UPCP, the dual of a TPCP map, then the 
Choi matrix enjoys the property 

TnC(r) = i, 

where Tr x denotes the partial trace over the first factor in C d ® C d . 

We shall compute the Jamiolkowski states of quasi-free TPCP maps Aa,b and V a ,b an d 
the Choi matrices of their duals. Both the Jamiolkowski state and the Choi matrix depend 
on the distinguished basis. A different choice of orthonormal basis in C d yields however a 
unitarily equivalent state and matrix. We can therefore equally well define J and C with 
respect to a naturally chosen basis in fermionic context. 

Jamiolkowski state of a TPCP map. 

Proposition 9. The Jamiolkowski state of a quasi-free TPCP map A a ,b or T a ,b is unitarily 
equivalent to a gauge-invariant quasi-free with symbol 

1/1 A \ 
2 \A* A* A + 2BJ 

for k A)B and 

If 1 -A \ 
2 \-A* A* A + 2B T J 

for T AjB . 

Proof. Consider first the case of a TPCP map Aa,b- 111 order to stay within the context of 
gauge-invariant quasi-free states we embed 21(5)) in the usual way in 21(53 © 53), identifying 
a(ip) with a(0 © if). The gauge-invariant quasi-free state on 21(53 © 53) with symbol 

, A A 

is pure and its marginal on 2l({0} © 53) is totally mixed, hence it can be used to construct 
the Jamiolkowski state. 

The algebra 21(53 © 53) can be decomposed as 21(53) © 21(53) but the factors cannot simply 
be chosen as 21(53 © 0) and 21(0 ©53), indeed, a(ipi © 0) anticommutes with a(0 © ip?). There 
exists in 21(53) an element such that 6* = 0, 2 = 1 and {0 , a(ip)} = 0. In fact, is up 
to a sign uniquely defined by these requirements and using the Fock space representation 
of the CAR, we easily see that = ±E(— 1). Let us now embed 21(53) in 21(53 © 53) as 
t(a(ip)) = ©i a(0 © ip) where ©i is the element in 21(53 © 0) just described. It is then easily 
checked that 21(53 © 53) decomposes into the tensor product of 21(53 © 0) and the embedded 
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algebra 21(55). Moreover, 

(id®A* AB ) (a((pi © ip 2 )) = (id®A* AB ) { a (Vi © 0) + a(0 © </? 2 )) 

= (id®A\ B ) (a((p! © 0) + 0i 0ia(O © ip 2 )) 
= (id <g)A^ B ) (a(y?i © 0) + 6i © i(a(p 2 ))) 
= (a(v?i©0) + ei©2(a(A^ 2 ))) 
= a((fi Q) Aip 2 ). 

In this way, we see that id®A A ,B is again quasi-free on 21(f) © $)) with defining operators 



A : = 



1 
A 



and B := 




Bj ■ 



(9) 



It is now obvious that the Jamiolkowski state will also be quasi-free with symbol 



CO 



i/ll 
2 > 1 1 



{ (id AA^ B ) (a*(y?i © ^M^i © ^2)) } 



M ! / 1 A 



Pi 

^2 



Remark that the positivity conditions for the symbol of the Jamiolkowski state precisely 
coincide with the positivity requirement < B < 1 — A* A for the TPCP map Aa,b- 

The computation for a map Ta,b is similar. It is now convenient to compose idAl^s 
with the local automorphism id A7 where 7(a(</?)) = a*(lp). In this way we remain within the 
class of gauge-invariant quasi-free states on the composite algebra 21(f) ©f)). The extended 
map is now of the form B with 



A 



1 

{A T y 



and 



B :-- 




B 



□ 



Choi matrix of a UPCP map. 



Proposition 10. The Choi matrix of a quasi-free UPCP map A* A B or T* A B is unitarily 
equivalent to 

'B- 1 - 1 B- L A* 
AB- 1 l + AB- 



Det (B) E 



I- 1 A*)] 



for A a ,b and 



Det (B) E 



f B- 1 -! B- l A J 
\{A J )*B- 1 1 + (A T )*B- 1 A T 



forT 



A,B ■ 
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Proof. On C d (g> C d the matrix £\ . (g> is equal to the dimension of the space times the 



projector on a maximally entangled vector. We should therefore compute 

2 d (id AA* a b )pq, 



where 



Q 



i i 
i i 



Using Proposition 6 we obtain 

2 d Det (1 - Q + (2Q - t)B) = Det 

and 

1 + A(1-Q + (2Q- t)B)-\2Q - 1)(A) 



1 -1 
2 \2B - 1 1 



Det (B) 



B- 1 - 1 B- X A* 
AB- 1 l + AB^A* 1 ' 



where A and B are as in (9). 

The computation for maps T* A B follows similar lines. Remark again that the positivity 
conditions for the Choi matrices coincide precisely with the conditions ensuring that 
and Ta b are CP. □ 
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